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ABSTRACT 

We investigate the range estimate between two wireless 
nodes without time stamps exchanging. Considering prac¬ 
tical aspects of oscillator clocks, we propose a new model 
for ranging in which the measurement errors include the 
sum of two distributions, namely, uniform and Gaussian. 
We then derive an approximate maximum likelihood es¬ 
timator (AMLE), which poses a difficult global optimiza¬ 
tion problem. To avoid the difficulty in solving the com¬ 
plex AMLE, we propose a simple estimator based on the 
method of moments. Numerical results show a promising 
performance for the proposed technique. 

1. INTRODUCTION 

Accurate distance estimates between two wireless nodes is 
a vital requirement for many applications, e.g., localizing 
the position of an unknown node |[Tl|2l . Among different 
approaches, the range estimate based on two-way time- 
of-arrival (TW-TOA) has attracted considerable attention 
in the literature and been proposed for practical applica¬ 
tions, e.g., in IEEE 802.15.4a lH. Traditionally, ranging 
based on TW-TOA involves exchanging the time stamps 
measured with respect to local clocks of different nodes. 
Namely, a master node initiates transmitting a signal at a 
certain time to a slave node and receives a response from 
the slave node after a delay corresponding to the distance 
between two nodes and a processing time, the so-called 
turn-around time a, at the slave node. The processing 
time can be estimated by the slave node, e.g., using a loop- 
back test HI, and then sent back to the master node. An¬ 
other approach to accurately estimate the range between 
two nodes is based on chirp-spread-spectrum that shows 
good performance in some situations iSIhl. 

It has been generally argued that ranging based on 
TW-TOA is severely affected by an imperfect clock 0 
Hill. An affine model, consisting of clock offset and 
clock skew parameters, is commonly considered to de¬ 
scribe the behavior of a clock of an oscillator HHIIl. Eor 
such a model, it is clearly seen that the distance estimate 
is mainly affected by an imperfect clock skew na. A 
number of studies tackled the positioning problem in the 
presence of unknown clock parameters in the past few 
years HiHia- 


In the literature, it is commonly assumed that the clock 
of an oscillator can be read continuously M131I16I . More¬ 
over, it is assumed that nodes are able to communicate 
with each other during the ranging period. In this study, 
we depart from these assumptions and consider sensor 
nodes as pure ranging devices. As it is common in prac¬ 
tice, a nominal value for the processing time is embedded 
in every node and once the slave node detects a signal, 
it responds according to the predefined turn-around time. 
Since in practice the clock of an oscillator can be read at 
certain times, e.g., at rising edges, we model the delay in 
detecting the arrival time by a uniform distribution. We, 
then, model the perturbation in ranging using the sum of 
two distributions, i.e., uniform and Gaussian, respectively, 
for modeling the delay in detecting the signal and the time- 
of-arrival estimation error. Eor such a model, we derive 
an approximate maximum likelihood estimator (AMLE), 
which poses a difficult global optimization problem due to 
nonconvexity of the AMLE objective function. We then 
propose a low complexity approach based on the method 
of moments (MOM) to estimate unknown parameters. The 
numerical results show a promising performance of the 
proposed approach, specially for high signal-to-noise ra¬ 
tios. In summary, the main contributions of this study are 
(a) a new ranging model for practical applications; (b) an 
approximate MLE (AMLE) for the proposed model; (c) 
a new low complexity technique based on the MOM for 
estimating the unknown parameters. 

2. SYSTEM MODEL 

Consider two master and slave nodes performing ranging 
using TW-TOA measurements. Eor the master node, we 
use a time-to-digital convertor (TDC) allowing us to mea¬ 
sure time stamp, which is not affected by imperfect clock 
of the master node. Note that the TDC is a passive de¬ 
vice; hence, it can not be used at the slave node to force 
the node to transmit the signal after certain processing de¬ 
lay. Eor details, see, e.g., El. In the literature, different 
models are used to describe a local clock. Among them, 
the popular one is an affine model expressed as ll^ fTbUTSll 

Cs{t) = eo + wt, ( 1 ) 

with Cs (t) as the local clock of the slave node with respect 
to the perfect time. In this model, Oq and w denote, respec- 


lively, the relative clock offset and the clock skew between 
the slave node and the reference time t. The model in ([T]i 
assumes that the clock of an oscillator can be read contin¬ 
uously, while in practice reading time happens at discrete 
instances, namely, at rising or falling edges. 

The relation between the clock skew and the oscillator 
frequency offset can be seen as follows. Let us consider 
the frequency fs of a slave node that deviates from nomi¬ 
nal frequency fo = I/Tq. The common imperfections are 
frequency offset and phase noise phenomena. Thus, we 
can model the frequency fs as 


Clearly, it is seen that the perturbation efc/2 + Uk has 
nonzero mean wTq /4, which depends on unknown clock 
skew. For large D, we may successfully neglect in (|6]l, 
especially for tiny To, and arrive at the traditional model 
of ranging considered in the literature. But, in general, 
the delay in detecting the arrival signal needs to be taken 
into account to have a more accurate model. In addition, 
a small value for the processing delay in slave node, small 
D, is preferable for some applications, e.g., for fast rang¬ 
ing. 

2.1. Maximum Likelihood Estimator 


fs = fo±Af + at), ( 2 ) 


where A/ denotes an offset and f{t) shows the perturba¬ 
tion, which is assumed to be zero-mean. Therefore, 


Ts = 


1 


fo±Af+m 


To{lTp) + m (3) 


where p A A///o is the deviation from ideal clock iw = 
1 =F p is called clock skew) and C,{t) = —f {t)/fo is known 
as clock jitter. Therefore, it is observed that the period Tg 
is a stochastic process with mean To(l =F p). In the rest 
of the paper, we assume that the jitter is small and can be 
neglected in the ranging process. 

We now investigate the ranging in the presence of 
clock imperfections at the slave node. In the fc-th round 
of performing TW-TOA, the master node sends a sig¬ 
nal to the slave node and receives a reply after a delay 
corresponding to the distance between the nodes and a 
processing delay T^, giving rise to the following model: 

d 

Zk = —I—I —k = 1,..., N (4) 
c 2 

where Uk is modeled by a zero-mean Gaussian random 
variable, i.e., Uk J\f{Q,cA) HKH], and c is the speed 
of propagation, d is the Euclidian distance between two 
nodes. We now model the delay as follows. In prac¬ 
tice, a nominal value for the delay is set into every node, 
say, = DTo, where D is an integer. It means that the 
slave node replies the detected signal after wDTo- There¬ 
fore, the total delay at slave node can be modeled by 

rf = wDTo + efe (5) 


We consider the vector of measurement z A \zi,..., zn^' 
and assume that and Uk are independent. In addition, it 
is assumed that nfc(or tk) and n^(or e^) are independent 
for k £. Thus, the probability density function (pdf) of 
the measurement vector 2 ; can be calculated as 


N ^wTo 

Pz{z;e) = n / PZki^klx, e)ps^{x)dx 



1 

y/^awTo 


exp 


/ {ak-x/2Y 
V 2a2 


dx (7) 


where ak A Zk — d/c — wDTo/2 and 6 A [d,w,a]. 
The integral in Q has no closed-form expression. For any 
X G K, we instead use the following approximation 1201 : 
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Using the above approximation, we obtain 
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where determines the delay in detecting the signal pres¬ 
ence. In fact, the arrived signal at the slave node may be 
detected after Ck, i.e., at the first rising edge of the clock. 
A natural way to model Cfc is to employ a uniform distri¬ 
bution Cfc ~ U{Q,wTo). In fact, for high-signal-to-noise 
ratios (SNRs), with high probability time-of-arrival detec¬ 
tion happens in the period that signal arrives. 

Replacing (|5]) into (01), we obtain the following model 
for ranging via TW-TOA: 


d wDTo Ck 

Zk = --\ -- V-zr+nk k = l,...,N. (6) 


where 13k = ak — wTq/2. Hence, we obtain the following 
approximation for the pdf of measurements in (I?]): 
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( 8 ) 



Fig. 1. AMLE cost function for |(a)| d and w for fixed a = 1/c, |(b)| d and tr for fixed w = 1.0001, and[^ w and a for 
fixed d = 50. 


We now obtain the AMLE expression by solving an opti¬ 
mization problem, expression ® shown at the top of the 
page. It is observed that the optimal estimator poses a dif¬ 
ficult global optimization problem. To get some feeling 
about the shape of the objective function, we fix one pa¬ 
rameter and plot the value of the objective function over 
two other parameters in Eig. [T] As it is observed, the ob¬ 
jective function for a fixed parameter has many maxima 
and it may also be discontinuous for some values. In the 
coming section, we propose a simple estimator based on 
the method of moments (MOM). 

2.2. A Low Complexity Estimator 


We first consider the relations between the unknown pa¬ 
rameters and the following statistics: 


Hi = Ezfe = - j i ] yjTf^ 


H2 = E(zfe - Hi) = cr + 


2 , (wTo)'^ 


48 


(9) 

( 10 ) 


Hi = - Hi) = 3cr -I- cr 


4 , 2(''"^o)^ , {wTo)"' 


1280 


( 11 ) 


The statistics pi, H 2 , Hi can be approximated by the 
means of ensemble averaging as follows: 


Ml 

M2 

Hi 


N 

N 


= S2 


= Si. 


where a = 1/1280 — 3/48^. Hence, an estimate of the 
clock skew is obtained as 


w = 





(13) 


We use the absolute value in (fOl l to prevent ambiguity due 
to noise in estimating the clock skew. 

Considering = (1 -I- p)^ ~ 1 + Kp, we can also 
estimate the clock skew from (fT^ as 


w = 1 + 


aT^ - Si + 35| 


(14) 


We can now estimate the distance and the variance cr^ as 




jwTo)^ 

48 


The mean of the estimator proposed in (fT4l) is 


Eu> = 1-1- 


aT4 


ES'4 -f 3ES'| 
4aT4 


(15) 

(16) 


(17) 


Using the law of large numbers, we have 


— ^ Ail - --f ( —-f-j wTo, N^qo. 


We now compute the expectations on the right-hand side 
of ( l20l i as follows. 


Erom ( [Tol l and (fTTT i. we obtain 

a{wTof‘ = Hi ~ 3/i2 


^*^4 = - ^4Si + Qzlsf - + St) 


( 12 ) 
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Fig. 2. Comparison between different approaches, |(a)| the 
mean of distance estimate and |(b)| the RMSE of the esti¬ 
mate. 


where Zk = Zk — /ii and Si = 1/N Since 

Zk are iid, we have 
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Fig. 3. The mean of the distance estimate versus the num¬ 
ber of samples N. 


3. SIMULATION RESULTS 

We compare the proposed technique with the traditional 
approach without clock skew compensation and a tech¬ 
nique based on approximating the turn-around time at the 
slave node via loop-back test and then sending back the 
estimate to the master node. In the simulation, we use 
d = 30 [m], /o = 100 MHz. We consider p = 0.0001, 
which corresponds to a frequency offset equal to 10 kHz. 
In the simulation, we set D = 10 and we collect TW-TOA 
measurements for 10 ms. To obtain the results, we run the 
algorithms for 1000 realizations of noise. 

In Fig. 12] we plot the mean of the distance estimate 
and the root-mean-square error (RMSE) of the estimate 
for different approaches. As it is observed, the proposed 
technique shows a considerable gain, especially for high 
SNRs. 

Fig 12] shows the mean of the estimate versus the num¬ 
ber of distance estimates, i.e., Eq. (fTSl l. N for fixed cr = 
0.1/c. As it is observed, after a sufficient number of sam¬ 
ples, e.g., 100 corresponding to 1 micro second, the esti¬ 
mate is very close to the theoretical value. 


4. CONCLUSIONS 


Likewise, we can show that 


ESl = {Ezlf+Oil/N). (19) 

Considering Ezf. = Str'^ -f ct^(wTo)^/ 8 -f (wTo)"‘/1280 
and Ezl = -\- (tcTo)^/48, it is verified that 


Ew = 1-1- 


1-f 


aT^ - a(wTo)^ 
4ar4 

aT^ - a{wTo)'^ 

4^1? 


+ 0{1/N) 

W, V > 1, 


( 20 ) 


In this paper, we have studied the distance estimate be¬ 
tween two nodes using TW-TOA. Considering practical 
aspects, we have modeled the range estimate using the 
sum of two distributions, uniform and Gaussian. We 
showed that the optimal estimator is complex and not 
easy to solve. Then, we proposed a low complexity tech¬ 
nique based on the method of moments. Numerical results 
show a promising performance for the proposed approach, 
especially for high signal-to-noise ratios. 
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